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ABSTRACT 



Aims. We announce the discovery of two new transiting planets, and provide their accurate initial characterization. 
Methods. First identified from the HATNet wide-field photometric survey, these candidate transiting planets were then followed-up 
with a variety of photometric observations. Determining the planetaiy nature of the objects and characterizing the parameters of the 
systems were mainly done with the SOPHIE spectrograph at the 1.93 m telescope at OHP and the TRES spectrograph at the 1.5 m 
telescope at FLWO. 

Results. HAT-P-42b and HAT-P-43b are typical hot Jupiters on circular orbits around early-G/late-F main sequence host stars, with 
periods of 4.641876±0.000032 and 3.332688±0.000016 days, masses of 0.975±0.126 and 0.660±0.083 Mj, and radii of 1.277±0.149 
and 1.283;^Q[j3^ Rj, respectively. These discoveries increase the sample of planets with measured mean densities, which is needed to 
constrain theories of planetary interiors and atmospheres. Moreover, their hosts are relatively bright (V < 13.5) facilitating further 
follow-up studies. 

Key words. Planetary systems - Techniques: photometry, radial velocity - Stars: individual: HAT-P-42, HAT-P-43 



L| 1 ■ Introduction physical properties of exoplanets, such as their mass-radius re- 

. !^ . lationships, atmospheric composition or orbital obliquity. They 

^ Currently more than 800 confirmed exoplanets are known, in- -^^ important first step toward a valuable comparison of 

^ c uding about 200 well-characterized transiting planets (lc. with ^^^^.^^^^^ These objects will be the targets of next space 

d planetary masses and radii measured to better than 10% accu- ^-^^-^^^ ^^^^ atmospheric characterization. 

■ ■ "racy). Transiting extrasolar planets are crucial for exploring the , ., , ^ , , ,• • 

RV surveys with well defined sample limits estimated an oc- 



* The photometric/spectroscopic data presented in this paper are currence rate of hot Jupiters in the Solai- neighborhood to be 
based in part on observations cartied out by the Hungarian-made Au- around one percent (Udry & Santos 2007; Wright et al. 2012). 
tomated Telescope Network, using telescopes operated at the Fred With a transit probability around ~ 10%, a huge number of stars 
Lawrence Whipple Observatory (FLWO) of the Smithsonian Astro- have to be monitored to detect such rare objects. This is why ma- 
physical Observatory (SAO), and at the SubmiUimeter Anray (SMA) joj-ity of transiting planets were discovered by photometric 
of SAO, by the Tillinghast Reflector 1.5 m telescope and the 1 .2 m tele- targeting tens of thousands of stai's per night. Ground- 
scope, both operated by SAO at FLWO, by the SOPHIE spectrograph , , J.UJ.. fU^ - J J 
I J 1 m . 1 . . ■ J TT / r> based surveys are biased toward the detection of short period and 
mounted on the 1.93 m telescope at Observatoire de Haute Provence, , 

France (rtins DDT-Dec201 1), by the Nordic Optical Telescope, oper- relatively large planets, i.e. hot Jupiters (with some exceptions, 
ated on the island of La Palma joinfly by Denmark, Finland, Iceland, e.g. HAT-P-1 1, Bakos et al. 2010). Complementanly, space- 
Norway, and Sweden, in the Spanish Observatorio del Roque de los based surveys (i.e. CoRoT and Kepler) have excelled in the de- 
Muchachos of the Instituto de Astrofisica de Canarias, and by the facil- tection of longer periods and smaller photometric transit depths, 
ities of the Las Cumbres Observatory Global Telescope. and hence smaller planet sizes. But, because the target stars are 
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Fig. 1. HATNet light curves of HAT-P-42 (top) and HAT-P-43 (bottom) phase folded with the transit period. In both cases we show two panels: 
the top shows the unbinned light curve, while the bottom shows the region zoomed-in on the transit, with dark filled circles for the light curve 
binned in phase with a binsize of 0.002. The solid line shows the model fit to the light curve. 



often relatively faint, radial-velocity follow-up to measure the 
planetary mass and confirm the planetary nature of the object is 
often not feasible for the space-based discoveries. 

We present here two new transiting hot Jupiters, first iden- 
tified from the Hungarian-made Automated Telescope Network 
wide-field photometric ground-based survey (HATNet; Bakos et 
al. 2004). Since 2006, HATNet has announced and pubHshed 41 
transiting exoplanets, i.e. ~ 20% of the well-characterized sam- 
ple. These two new typical hot Jupiters, HAT-P-42b and HAT-P- 
43b, were partially confirmed thanks to spectroscopic observa- 
tions made with SOPHIE mounted at the 1 .93 m telescope at the 
Observatoire de Haute-Provence (OHP) in France and TRES at 
the 1 .5 m telescope at FLWO. 

The layout of the paper is as follows. Section 2 presents the 
diverse photometric and spectroscopic observations that lead to 
the detection and characterization of these systems. In the third 
section, we derive the stellar parameters and the planetary orbits. 
The characteristics of these two hot Jupiters are discussed in the 
last section. 



2. Observations 

2.1. Photometric Detection 

HAT-P-42 and HAT-P-43 were identified as candidate transit- 
ing planets based on photometric observations conducted by the 
HATNet survey (Bakos et al. 2004). These observations are 
summarized, together with the follow-up photometric observa- 
tions, in Table[T] The data were processed and searched for tran- 
sits following the procedure of Bakos et al. (2010), see also 
Kovacs et al. (2002, 2005). Figure [U shows the folded HATNet 
light curves for both systems. We give cross-identifications, and 
catalog photometry on an absolute scale for each system in Ta- 
ble [T] 

2.2. Follow-up 

We conducted follow-up spectroscopic and photometric obser- 
vations of both systems to confirm their planetary natures and 
determine their physical properties. See Latham et al. (2009) for 
a more detailed description of our follow-up procedure. Below 
we provide specific details relevant to the discoveries of HAT-P- 
42b and HAT-P-43b. 



Table 2. High-precision differential photometry of HAT-P-42. 



BJD" 


Mag" 


C Mag 


Mag(oiig)'^ Filter 


(2,400,000h-) 






55882.87909 


0.01318 


0.00170 


11.46340 / 


55882.88049 


0.01025 


0.00191 


11.46150 / 


55882.88128 


0.01064 


0.00192 


11.45990 i 


55882.88208 


0.00410 


0.00192 


11.45190 i 


55882.88303 


0.00898 


0.00173 


11.45770 i 


55882.88389 


0.01010 


0.00180 


11.46170 i 


55882.88473 


0.00376 


0.00187 


11.45210 / 


55882.88558 


0.00885 


0.00187 


11.45610 / 


55882.88644 


0.01057 


0.00188 


11.45900 i 


55882.88729 


0.00559 


0.00168 


11.45500 i 



Notes. This table is available in a machine-readable form in the on- 
line journal. A portion is shown here for guidance regarding its form 
and content. '^'^ Barycentric Julian Date calculated directly from UTC, 
without correction for leap seconds. The out-of-transit level has been 
subtracted. These magnitudes have been subjected to the EPD and TFA 
procedures, carried out simultaneously with the transit fit. Raw mag- 
nitude values without application of the EPD or TFA procedures. 

2.2.1. Reconnaissance Spectroscopy 

We obtained initial high-resolution, low S/N "reconnaissance" 
spectra of HAT-P-42 and HAT-P-43 using the Tillinghast Re- 
flector Echelle Spectrograph (TRES; Furesz 2008) on the 1.5 m 
Tillinghast Reflector at the Fred Lawrence Whipple Observatory 
(FLWO) in AZ, USA, and using the Fibre-fed Echelle Spec- 
trograph (FIES) on the 2.5 m Nordic Optical Telescope (NOT) 
on the island of La Palma, Spain (Djupvik & Andersen 2010). 
These observations are summarized in Table 2] The data were 
analyzed via cross-correlation against synthetic templates as de- 
scribed by Buchhave et al. (2010) and Quinn et al. (2012). Both 
stars were found to be slowly rotating G dwarfs, with no evi- 
dence for composite spectra, and with low velocity variations 
indicating that neither object is an undiluted ecUpsing binary star 
system. 

2.2.2. Photometry 

Photometric follow-up observations were obtained with Kepler- 
cam on the FLWO 1 .2 m, the Spectral CCD on the 2.0 m FauUces 
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Fig. 2. Follow-up light curves for HAT-P-42 (left) and HAT-P-43 (right). For HAT-P-42 all observations were obtained with Keplercam on the 
FLWO 1.2 m telescope. For HAT-P-43 observations were obtained with Keplercam on the FLWO 1.2 m (indicated by FLWO in the figure), the 
Spectral CCD on the FTS 2.0 m (FTS), and CCD imager on the BOS 0.8 m (BOS). The Keplercam light curves have been corrected for trends 
during the modeling. The BOS and FTS light curves have also been corrected for trends, though for these light curves we only apply the External 
Parameter Decorrelation (EPD) procedure, and not the Trend Filtering Algorithm (TFA) (see for details Bakos et al. 2010). The dates of the events 
are indicated. Curves after the first are displaced vertically for clarity. Our best fit from the global modeling is shown by the solid lines. Residuals 
from the fits are displayed at the bottom, in the same order as the top curves. The error bars represent the photon and background shot noise, plus 
the readout noise. They are only plotted on the residuals for reason of readability. 



Table 1. Summary of Photometric Observations. 



Instrument/Field'' 


Date(s) 


Number of Images 


Approx. Cadence (min) 


Filter 


HAT-P-42 










HAT-6/G366 


2OION0V-2OII Apr 


5325 


3.7 


Sloan r 


HAT-9/G366 


2OION0V-2OII Apr 


4859 


3.7 


Sloan r 


Keplercam 


2011 Nov 16 


176 


1.2 


Sloan i 


Keplercam 


2012 Jan 6 


310 


1.0 


Sloan i 


Keplercam 


2012 Feb 3 


93 


1.2 


Sloan i 


Keplercam 


2012 Mar 2 


278 


1.1 


Sloan i 


HAT-P-43 










HAT-5/G317 


2010 Nov-2011 Apr 


4174 


3.7 


Sloan r 


HAT-8/G317 


2OION0V-2OII Apr 


4300 


3.7 


Sloan r 


HAT-7/G365 


2010 Nov-2011 May 


8262 


3.7 


Sloan r 


HAT-8/G365 


2011 Apr-2011 May 


464 


3.7 


Sloan r 


Keplercam 


2011 Nov 17 


68 


3.2 


Sloan i 


Keplercam 


2012 Jan 13 


87 


2.2 


Sloan i 


FTS 


2012-Feb-16 


226 


1.0 


Sloan i 


Keplercam 


2012 Mar 3 


141 


2.2 


Sloan i 


BOS 


2012 Mar 3 


109 


2.4 


Sloan r 


Keplercam 


2012 Mar 23 


125 


2.2 


Sloan i 



Notes. HAT-5, -6, -7, and -10 are located at FLWO in AZ, USA. HAT-8 and -9 are located at Mauna Kea Observatory in HI, USA. 
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Fig. 3. High precision RV measurements and bisector spans (BS) for HAT-P-42 (left) and HAT-P-43 (rigiit). Filled circles show 

OHP 1.93m/SOPHIE observations, while open triangles show FLWO 1.5m/TRES observations. The top panels show the RV measurements 
as a function of orbital phase, along with our best-fit circular orbit model. Zero phase corresponds to the time of mid-transit. The center-of-mass 
velocity has been subtracted for each system. The middle panels show velocity 0-C residuals from the best fit circular orbit model. The bottom 
panels show the BS values, with the mean values subtracted. Note the different vertical scales of the panels. 
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Fig. 4. Bisector spans (BIS) as a function of the high precision RV measurements for HAT-P-42 (left) and HAT-P-43 (right) from the 
OHP 1.93 m/SOPHIE observations. In both panels the scale is the same in the x and y axes. The mean values of the BS have been subtracted for 
each system. 



Telescope South (FTS) at Siding Springs Observatory in Aus- 
tralia, and with the CCD imager on the 0.8 m Byrne Observatory 
at Sedgwick (BOS) telescope at Sedgwick Reserve in the Santa 
Ynez Valley, CA, USA. These observations are summarized in 
Table [1] Both FTS and BOS are operated by the Las Cumbres 
Observatory Global Telescope (LCOGlQ; Shporer et al. 2011; 



' http://lcogt.net 



Brown et al. 2013, in preparation). The Keplercam observations 
were reduced to light curves following the procedure of Bakos 
et al. (2010), while the FTS and BOS observations were reduced 
following the procedure of Fulton et al. (201 1). 

Figure |2] shows the photometric follow-up light curves for 
both systems together with our best-fit models. HAT-P-42 and 
HAT-P-43 have photometric transit depths of 7.6 and 16.5 mmag. 
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Fig. 5. Model isochrones from ? for the metallicities of HAT-P-42 (left) and HAT-P-43 (right) and ages of 1 to 13 Gyr in 1 Gyr increments (left 
to right). The adopted values of Teff* and a/^* are shown together with their la and Icr confidence ellipsoids. In each plot the initial values of 
Teff* and a/R^, from the first SPC and light curve analyses are represented with a triangle. 
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Fig. 6. Mass radius plot of transiting extrasolar planets with determined mass. HAT-P-42 and HAT-P-43 are marked with the large filled circles. 



respectively. The data are provided in machine-readable format 
in Table |2]for HAT-P-42, and in Table[3]for HAT-P-43. 



2.2.3. Confirmation Spectroscopy 

We obtained high-resolution, high-S/N spectra of HAT-P-42 and 
HAT-P-43 using the SOPHIE spectrograph on the 1.93 m tele- 
scope at OHP (Bouchy et al. 2009), and of HAT-P-42 using 
TRES on the FLWO 1.5 m telescope. These observations are 
summarized in Table |4] The TRES observations of HAT-P-42 
were simply a continuation of the first two "reconnaissance" 
TRES observations of this target which were already of sufficient 
S/N to detect the ~ 100 ms^' orbital variation. Of the 13 TRES 
measurements, 4 with too low S/N (i.e. precision > 100ms"') 
were excluded from the orbital analysis. While we have previ- 
ously used SOPHIE to confirm HATNet planets (e.g. Bakos et al. 
2007, Shporer et al. 2009), there have been significant changes 
to the instrument and reduction procedure. These significantly 
increased the RV accuracy (Pemichot et al. 2011, Bouchy et 



al. 2012), thus we provide a description of our observing and 
reduction procedures below. 

We observed HAT-P-42 and HAT-P-43 with the High Effi- 
ciency fibers, yielding a resolving power of R-A/AA » 39 000 
(at 550 nm). The spectra were correlated with a G2 numerical 
mask in order to calculate the cross-correlation function (CCF). 
This CCF is fitted by a Gaussian and its center gives the RV 
measurement (Baranne et al. 1996, Pepe et al. 2002). The bi- 
sector spans values were calculated as described in Queioz et al. 
(2001). The error bars were derived from the CCF, following 
the procedure in Boisse et al. (2010), yielding a mean precision 
of ~ 6ms-' for HAT-P-42 and ~ 10m s ' for HAT-P-43. Since 
most of our spectra were polluted by Moon light, we corrected 
them following the method presented in Hebrard et al. (2008), 
and increased the error bars as determined empirically, leading 
to a mean precision of ~ 40ms-' for HAT-P-42 and ~ 30ms-' 
for HAT-P-43. A thorium-argon calibration was done each two 
hours to monitor the instrumental drift mainly due to small vari- 
ation of temperature and pressure in the spectrograph tank. SO- 
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Fig. 7. Radius of transiting extrasolar planets as a function of their equilibrium temperature (T^,). HAT-P-42 and HAT-P-43 are marked with big 
filled circles. Their respective T^q values are calculated by assuming zero Bond albedos, and also that the flux is re-radiated from the full planet 
surface. It is not the case for all planets and we advise the reader that the T^^ sample is not homogeneous. 

Table 4. Summary of Spectroscopic Observations. 



Instrument 


Date(s) 


Number of 


Resolution 


Wavelength 


Velocity" 






Observations 


Ai ^A 


Coverage [A] 


Precision 


HAT-P-42 












FLWO 1.5m/TRES 


2011 Oct-2012Jan 


13 


44,000 


3900-8900 


~ (30 + 35)'' ms-i 


OHP 1.93m/SOPHIE 


2011 Dec 7-13 


5 


39,000 


3900-6900 


~ (6 + 35f ms-i 


HAT-P-43 












NOT 2.5 m/FIES 


2011 Oct 24-27 


2 


46,000 


3700-7300 


~ lOOms-i 


FLWO 1.5m/TRES 


2011 Nov-2012Jan 


2 


44,000 


3900-8900 


~ 100ms"' 


OHP 1.93m/SOPHIE 


2011 Dec 6-13 


6 


39,000 


3900-6900 


~ (10 + 25)'= ms-i 



Notes. We give the median velocity precision of the observations, and not the best precision that can be attained with the given instrument. 
The FIES and TRES observations of HAT-P-43 were done in reconnaissance mode (lower S/N) which is why the velocities are less precise than 
the SOPHIE observations. A precision of 30ms"', plus 35ms"' added in quadrature to achieve a. x' per degree of freedom of unity for the 
Keplerian fit of HAT-P-42. Values in the parentheses include the sum of the theoretical noise due to the instrument (6 and 10m s"', respectively), 
plus empirical errors taking into account the correction due to Moon light in the spectra. 



PHIE is environmentally stabilized and this intrinsic drift is less 
than 2m s"' per hour Each measurement was also corrected for 
this drift. 

We computed the CCF using masks corresponding to differ- 
ent spectral types (FO, G2, K5, M5) and found no significant 
variation in the RV semi-amplitude. If a planet induces a RV 
shift, the amplitude of the variation should remain constant re- 
gardless of the spectral lines used to measure the RV. In contrast, 
for the case of a blend consisting of a bright target and a faint 
binary having different spectral types, the relative contributions 
to the CCF from the different components will change when a 
different mask is used (Santos et al. 2002, Collier Cameron et al. 
2007). This result thus favors the planetary hypothesis. 

Figure [3] shows the high-precision RV curves for both sys- 
tems, together with our best-fit circular-orbit models. Figure |4] 
shows the bisector spans as a function of the RV from the SO- 
PHIE spectra. The data are provided in Table |5] for HAT-P-42, 
and in Table|6]for HAT-P-43. 



3. Analysis 

To rule out the possibility that either object might actually be a 
blended stellar eclipsing binary system, we conducted a blend 
analysis similar to that done in Hartman et al. (2012). We find 
that in both cases we cannot rule out blend scenarios based on 
the photometric observations alone, however all scenarios which 
fit the photometric data predict significant RV and BS variations 
(greater than 300m s"' for HAT-P-42 and greater than 1 km s"' 
for HAT-P-43). Such variations are ruled out by our spectro- 
scopic observations. We conclude that both objects are transiting 
planet systems. 

Having confirmed the planetary nature of both systems, we 
proceeded with their analysis following the methods described 
in Bakos et al. (2010), with some modifications as described 
in Hartman et al. (2012). Briefly, this consists of: (1) in- 
fening the stellar atmospheric parameters from the available 
high-resolution spectra (we used the TRES spectra of HAT-P- 
42, and the TRES and FIES spectra of HAT-P-43, together with 
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Table 3. High-precision differential photometry of HAT-P-43. 



HAT-P-42b and HAT-P-43b 

Table 6. Relative radial velocities and bisector spans of HAT-P-43. 



BJD" 






Mag" 




O" Mag 


Mag(orig)'^ Filter 


(2,400,0004-) 










55883 


89634 





00397 





00167 


12.40490 / 


55883 


89856 





00002 





00157 


12.40130 / 


55883 


90082 


-0 


00162 





00153 


12.39990 / 


55883 


90306 





00397 





00150 


12.40600 i 


55883 


90530 


-0 


00219 





00158 


12.39980 i 


55883 


90753 





00141 





00150 


12.40250 / 


55883 


90976 


-0 


00085 





00155 


12.40180 / 


55883 


91201 


-0 


00087 





00153 


12.40340 / 


55883 


91426 





00004 





00151 


12.40220 / 


55883 


91650 


-0 


00165 





00166 


12.39820 / 



Notes. This table is available in a machine-readable form in the on- 
line journal. A portion is shown here for guidance regarding its form 
and content. Barycentric Julian Date calculated directly from UTC, 
without correction for leap seconds. The out-of-transit level has been 
subtracted. These magnitudes have been subjected to the EPD and, in 
some cases, TFA procedures, carried out simultaneously with the tran- 
sit fit. The TFA procedure has not been applied to observations with 
2455973 < BJD < 2455975 or to observations in the r filter with 
2455990 < BJD < 2455991. Raw magnitude values without ap- 
plication of the EPD or TFA procedures. 

Table 5. Relative radial velocities and bisector spans of HAT-P-42. 



BJD^ 


RV" 




Els'* 


Phase 


(2,454,000-H) 


(ms-i) 


(ms-i) 


(ms-i) 




1889.03539 


-47.4 


28.4 




0.322 


1900.05380 


159.1 


33.8 




0.696 


1901.01864 


27.1 


32.4 




0.904 


1902.02682 


-105.1 


34.6 




0.121 


1902.71983 


-112.1 


13.0 


-43.0 


0.270 


1902.93220 


-75.4 


25.4 




0.316 


1903.67805 


8.9 


31.0 


-27.0 


0.477 


1904.03850 


-14.2 


25.4 




0.554 


1904.66696 


97.9 


38.0 


5.0 


0.690 


1905.02718 


75.3 


29.9 




0.767 


1906.60475 


-46.1 


44.0 


8.0 


0.107 


1907.02711 


-33.4 


25.9 




0.198 


1908.63608 


71.9 


69.0 


-35.0 


0.545 


1939.89586 


-166.5 


33.4 




0.279 



Notes. Observations without a BS value were obtained with the TRES 
instrument, observations with a BS value were obtained with SOPHIE. 

Barycentric Julian Date calculated directly from UTC, without cor- 
rection for leap seconds. A zero-level y velocity, fitted independently 
for each instrument, has been subtracted from these measurements. For 
reference, the SOPHIE y velocity from the fit was 20.561±0.014km s"', 
which is an estimate of the true recession velocity of the center of mass 
of the system. For TRES the value is 0.014±0.015km s"' . As an artifact 
of the reduction procedure the TRES RVs were measured relative to an 
arbitrary zero-point, the fitted y velocity in this case has no physical 
meaning. Internal errors excluding any component of astrophysical 
jitter. We use the relation ctbs = Zctrv to estimate the BS uncertain- 
ties. 



the Stellar Parameter Classification (SPC) method of Buchhave 
et al. (2012); (2) conducting a global Markov-Chain Monte 
Carlo (MCMC)-based modeling of the available photometric 
light curves and RVs (we fix the limb darkening coefficients us- 
ing the tables by Claret 2004); (3) using the spectroscopically 
inferred stellar effective temperatures and metallicities, together 



BJD^ 


RV" 


ctrv" 


BIS" 


Phase 


(2,454,000+) 


(ms-i) 


(ms-i) 


(ms-i) 




1901.68768 


-85.8 


15.0 


-36.0 


0.289 


1902.67245 


40.3 


15.0 


-77.0 


0.585 


1903.62875 


42.3 


31.0 


12.0 


0.872 


1904.62601 


-94.8 


27.0 


-18.0 


0.171 


1906.57509 


111.3 


28.0 


-26.0 


0.756 


1908.60322 


-28.8 


31.0 


35.0 


0.364 



Notes. All observations listed were obtained with SOPHIE. Barycen- 
tric Julian Date calculated directly from UTC, without correction for 
leap seconds. For reference, the SOPHIE y velocity from the fit was 
-5.082 ± 0.009km s"'. Internal errors excluding any component of 
astrophysical jitter. We use the relation trgs = 2crRv to estimate the 
BS uncertainties. 

with the stellar densities determined from the light curve mod- 
eling, and the Yonsei-Yale theoretical stellar evolution models 
Yi et al. (2001), to determine the stellar masses, radii and ages, 
as well as the planetary parameters (e.g. mass and radius) which 
depend on these values (Figure |5]l; (4) re-analyzing the high- 
resolution spectra fixing the stellar surface gravities to the values 
found in (3), and then re-iterating steps (2) and (3). 

In conducting the analysis we inflated the TRES RV errors 
for HAT-P-42 by adding in quadrature a value of 35 m s ' to the 
internal errors. This "jitter" is needed to achieve a per degree 
of freedom of unity for the TRES RVs in the best-fit model. We 
do not add jitter to the SOPHIE RVs of HAT-P-42 or HAT-P-43 
because in both cases per degree of freedom is less than one 
for these RVs, indicating that the SOPHIE formal uncertainties 
may be overestimated. 

For both systems we conducted the analysis fixing the eccen- 
tricities to zero, as well as allowing the eccentricities to vary. For 
each system we find that the eccentricity is consistent with zero 
(the 95% confidence upper limits on the eccentricity are e < 0.2 
for HAT-P-42b, and e < 0.29 for HAT-P-43b). Following the 
suggestion of Anderson et al. (2012) we adopt the parameter val- 
ues associated with the fixed circular orbits. The adopted stellar 
parameters are given in Table Q while the adopted planetary pa- 
rameters are given in Table |9] For completeness we also provide 
the parameters which result when the eccentricities are allowed 
to be non-zero in Tables [8] and [10] Note that all the eccentric 
parameters are within 1 cr of the circular orbit values. 

We searched for sinusoidal signals and additional transits in 
the light curves. Neither present a significant peak in the L-S 
periodogram (false alarm probability of 50% or more). After 
removing the detected in-transit data, neither light curves shows 
a transit signal that would had passed our automated or by-eye 
selections. 

4. Discussion 

We have presented the discovery, confirmation and characteriza- 
tion of two new transiting planets. 

HAT-P-42b and HAT-P-43b are inflated hot Jupiters with 
P = 4.641876 + 0.000032 and 3.332688 ± 0.000016 days, 
Mp = 0.975+0.126 and0.660+0.083Mj,and7?p = 1.277+0.149 
and 1.283^0 034 ^J' respectively. Assuming zero albedo and full 
heat redistribution, and using the stellar and planetary parame- 
ters determined by our analysis, their equilibrium temperatures 
are in the same regime, i.e. T^q = 1427 + 58 and 1361 + 24 K, re- 
spectively. HAT-P-42b has larger uncertainties on its deiived pa- 
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Table 7. Adopted stellar parameters for HAT-P-42 and HAT-P-43 assuming circular orbits. 



HAT-P-42 HAT-P-43 

Parameter Value Value Source 

Astrometric properties 

GSCID GSC 0232-01451 GSC 0801-00608 

2MASS-ID .... 2MASS 09012265+0605500 2MASS 08354217+1012239 

R.A. (J2000) . . . 09'>01"22.66s 08h35"42.18s 2MASS 

Dec. (J2000) . . . +06°05'50.0" +10°12'24.0" 2MASS 

jUra (masyr-i) . -6.2±1.9 -10.3 + 2.6 UCAC4 

//Dec (masyr-i). -29.3 + 2.0 -16.0 + 3.2 UCAC4 

Spectroscopic properties 

reff*(K) 5743 ±50 5645 + 74 SVC 

[Fe/H] 0.27 + 0.08 0.23 + 0.08 SPC 

wsin/(kms-i).. 3.5 + 0.5 2.4 + 0.5 SPC 

rRv (kms-i) ... 20.26 + 0.15'' -4.86 + 0.15'' TRES 
Photometric properties 

B(mag) 12.827 + 0.020 14.120 + 0.060 APASS 

y(mag) 12.168 + 0.030 13.356 + 0.030 APASS 

7 (mag) 10.960 + 0.023 12.146 + 0.021 2MASS 

//(mag) 10.677 + 0.024 11.809 + 0.029 2MASS 

(mag) 10.626 + 0.026 11.764 + 0.023 2MASS 

Derived properties 

M^(Ms) 1.179 + 0.067 l-048+°°^^ YY+a/R^+SPC^ 

Ri,(RQ) 1.528 + 0.136 l-104!°-°^^ YY+a/R^+SPC 

log0*(cgs).... 4.14 + 0.07 4.37 + 0.02 YY+a/R^+SPC 

L^(Lq) 2.27 + 0.42 1.12 + 0.10 YY+a/R^+SPC 

My (mag) 3.94 + 0.21 4.73 + 0.11 YY+a/R^+SPC 

M^(mag,ESO). 2.38 + 0.20 3.11 + 0.06 YY+a/R^+SPC 

Age(Gyr) 5.1+i-^ 5.7+|-^ YY+a//;*+SPC 

Av (msigf 0.000 + 0.020 0.000 + 0.019 YY+a/R^+SPC 

Distance (pc) . . . 447 + 40 543+ 18 YY+a/R^+SPC 



Notes. SPC = "Stellar Parameter Classification" method, described by Buchhave et al. (2012), which derives stellar atmospheric parameters 
from high-resolution spectra. These parameters rely primarily on SPC, but have a small dependence also on the iterative analysis incorporating the 
isochrone search and global modeling of the data, as described in the text. These velocities corresponds to an absolute scale references thanks 
to nightly observations of HD 182488 and correction of gravitational redshift of the Sun. *' YY+a/i?*+SPC = Based on the YY isochrones (Yi et 
al. 2001), a/R-), as a luminosity indicator, and the SPC results. V band extinction determined by comparing the measured 2MASS and APASS 
photometry for the star to the expected magnitudes from the YY+a/R^,+SPC model for the star. We use the Cardelli et al. (1989) extinction law. 

Table 8. Derived stellar parameters for HAT-P-42 and HAT-P-43 allowing nonzero eccentricty. 





HAT-P-42 


HAT-P-43 


Parameter 


Value 


Value 


(Mo) 


1.185 + 0.070 


1.044 + 0.050 


R*(Re) 


1.497 + 0.190 


1 151+°-2«' 


Iog0* (cgs) .... 


4.16 + 0.09 


4.33 + 0.09 


(Le) 


7 70+0-68 
^■^^-0.48 


1 1Q+0.50 
^•^^-0.18 


My (mag) 


3.97 + 0.28 


4.66 + 0.27 


Mk (mag,ESO) . 


2.42 + 0.28 


3.02 + 0.26 


Age (Gyr) 




6.5+2-^ 


Av (magf 


0.000 + a023 


0.000 + aoi4 


Distance (pc) . . . 


439 + 56 


563!!? 



Notes. Quantities and abbreviations are as in Table [7] which gives our adopted values, determined assuming circular orbits. We do not list 
parameters that are independent of the eccentricity. 



rameters as the transit depth is shallower due to the larger radius 
of the slightly more evolved host. In Fig. |6] and Fig.|7]we ploQ 
these planets in the mass-radius and Teq-radius diagram, show- 
ing that their positions are typical compared to the sample of 

^ We used the up-to-date catalogue of the transiting planets TEPCat 
|http : //www ■ astro . keele . ac . uk/ jkt/tepcat/. 



well-characterized transiting exoplanets. Within 1 cr, HAT-P-42b 
and HAT-P-43b have the same radius and Teq. They mainly dif- 
fer by their periods (semi-major axes) and by their masses (den- 
sities). HAT-P-43b is very similar to HAT-P-4b (M,, = 0.672Mj, 
P = 3.05d, Rp = 1. 27 R3, Kovacs et al. 2007) with HAT-P- 
4 being a little more massive, hotter and more evolved star. 
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Boisse et al.: HAT-P-42b and HAT-P-43b 
Table 9. Adopted orbital and planetary parameters for HAT-P-42b and HAT-P-43b assuming circular orbits. 



HAT-P-42 HAT-P-43 

Parameter Value Value 
Light curve parameters 

P(days) 4.641876 + 0.000032 3.332688 ± 0.000016 

T, (BJDT 2455952.52606 + 0.00077 2455997.37105 ± 0.00033 

ri4(daysf 0.1680 + 0.0038 0.1355 + 0.0011 

Ti2 = Tu (days)" .... 0.0193 + 0.0036 0.0150 + 0.0008 

a/R. 8.09+0™ 8 63+0.12 

^/R^^ 13.40 + 0.14 16.60 + 0.09 

Rp/R^ 0.0858 + 0.0032 0.1194 + 0.0018 

^2 Q 332+0.084 039+" "^^ 

b = acosi/R^ 0.577!°o^^ 0.197+|;';f,° 

/ (deg) 85.9!';^ 88.7 + 0.7 

Limb-darkening coefficients'^ 

ci,/ (linear term) .... 0.2863 0.3003 

C2, / (quadratic term) . 0.3274 0.3180 

RV parameters 

K(ms-^) 105.8+ 13.0 87.5+ 10.8 

e (fixed) (fixed) 

Planetary parameters 

MpiMi) 0.975 + 0.126 0.660 + 0.083 

Rp{Ri) 1.277 + 0.149 1-283+;};;^^ 

C{Mp,Rpf 0.18 0.07' 

Pp {gcm.-^) 0-58!oi5 0.38 + 0.06 

logg^, (cgs) 3.17 + 0.11 2.99 + 0.06 

a(AU) 0.0575 + 0.0011 0.0443!°°™4 

req(K)'= 1427 + 58 1361 + 24 

0f 0.074+o°|^ 0.043 + 0.006 

<F)(10''ergs-'cm-^)g 9.36 + 1.52 7.75 + 0.54 



Notes. Reported times are in Barycentric Julian Date calculated directly from UTC, without correction for leap seconds. Tc'. Reference epoch 
of mid transit that minimizes the correlation with the orbital period. T^: total transit duration, time between first to last contact; T[2 = T^^: 
ingress/egress time, time between first and second, or third and fourth contact. Reciprocal of the half duration of the transit used as a jump 
parameter in our MCMC analysis in place of ajR-i,. It is related to ajR-,, by the expression (IR-,, = a/Rj,{2jT{l + esin(jj))/(PVl - Vl - e^) 
Bakos et al. (2010). Values for a quadratic law, adopted from the tabulations by Claret (2004) according to the spectroscopic (SPC) parameters 
listed in Table |7] Correlation coefficient between the planetary mass M,, and radius Rp. Planet equilibrium temperature averaged over the 
orbit, calculated assuming a Bond albedo of zero, and that flux is re-radiated from the full planet surface. The Safronov number is given by 
= T(Vesc/Vorb)' = (alR„){M„IM^) (see Hansen & B arman 2007). Incoming flux per unit surface area, averaged over the orbit. 



with similar metallicity ([Fe/H] = 0.24), and to OGLE-Tr- 
10b {Mp = 0.63Mj, P = 3.1d, = 1.25/?j, Torres et al. 
2008). HAT-P-42b is resembling to CoRoT- 1 9b (Mp = 1.107Mj, 
P = 3.9d, Rp = 1.29/?j, Guentheret al. 2012). HAT-P-13b is 
the other closest analogous to HAT-P-42b and HAT-P-43b but 
in the intermediate mass-regime {Mp - 0.85 IMj, P - 2.9d, 
Rp = 1.287?:, Winn et al. 2010, Bakos et al. 2009). The stellar 
hosts of all the previously quoted planets are slow rotators, with 
late F/early G spectral type, and with enhanced metallicity with 
respect to the Sun (except for CoRoT-19, which is solar). 

HAT-P-42b and HAT-P-43b are inflated hot Jupiters com- 
pared to models of coreless giant planets. The models of Fort- 
ney et al. (2007) predict a maximal radius of ~ 1 . \Rj for both 
planet (see their Table 4 and Figure 7) considering the system 
age, planet distance and mass (maximal in the meaning of a core 
free planet of pure H-He). Laughlin et al. (201 1) derived a sim- 
ple fitting relation from the Bodenheimer et al. (2003) models 
to infer the expected radius of a H-He composition planets for 
different masses and insolation. From the Laughlin et al. (201 1) 
relations, we derived an predicted radius of 1.19 and \.\lRj, re- 
spectively for HAT-P-42b and HAT-P-43b, still at the 1-cr and 
3-0" level from our observations. 



Enoch et al. (2012) determined an empirical relation be- 
tween Rp, Teq, and the semi-major axis a (see their Eq. 9). This 
leads to an estimated radius of 1.325 + 0.054 /?j for HAT-P-42b, 
matching the observed radius within the error bars. The same 
relations predict a radius of l.im^^Qf^^Rj for HAT-P-43b, which 
is 3-cr below our measured value. Considering also the results 
observed for the highly inflated hot Jupiters presented by Hart- 
man et al. (2012), it is possible that other parameters than Teq 
and a play a role in the planetary radius in the mass domain 
where the equation was derived for 0.5 < Mp < 2Mj. Indepen- 
dently, Beky et al. (2011) derived a relation to determine the 
radius for the planets with 0.3 < Mp < 0.8Mj that we could 
apply to HAT-P-43b. Their equation leads to an inferred ra- 
dius of 1.133 + 0.044 As was already remarked for HAT- 
P-39b/40b/41b by Hartman et al. (2012), the predicted radius is 
smaller than the measured value. 

Another caveat in using the above relations is the assumption 
of zero Bond albedo and full redistribution of the heat when cal- 
culating Teq, whereas in reality these assumptions are unlikely to 
be true for all of the planets. 

Few measurements per orbit lead to poor constraints on a 
small eccentricity. We followed the recommendation of Ander- 
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Table 10. Orbital and planetary parameters for HAT-P-42b and HAT-P-43b allowing eccentric orbits. 



HAT-P-42 HAT-P-43 

Parameter Value Value 

Light curve parameters 

o/R* 8.28^Qg3 8.25^jq| 

^/R* 13.41 ±0.14 16.62 ±0.09 

!(deg) 86.2 ± 1.2 88.5+','-[ 

RV parameters 

K(ms-^) 104.4 ± 13.8 89.7 ±13.0 

Vecosw -0.060±0.158 -0.123±0.181 

Vesinw -0.055 ± 0.226 0.146 ±0.228 

ecosw -0.012±0.052 -0.031!°°^5 

esinw -0.008 ±0.084 0.035+[; i^° 

e 0.067 ±0.062 0.091 ±0.094 

w(deg) 208 ±91 135 ± 80 

Secondary eclipse parameters 

(BJD) 2455954.812 ±0.155 2455995.639 ± 0.161 

T,i4 0.1660 ±0.0176 0.1451 ± 0.0341 

r,,i2 0.0179 ± 0.0090 0.0164 ± 0.0056 

Planetary parameters 

MpiMj) 0.957 ±0.137 0.674 ±0.097 

Rp(Ri) 1.250 ±0.185 ^-^"^Kfrn 

CiMp,Rp) 0.33 0.50 

Ppigcm-^) 0.61+5;]^ 0.34±0.10 

loggpicgs) 3.18 ±0.12 2.96!°?^ 

a (AU) 0.0576 ± 0.001 1 0.0443 ± 0.0007 

req(K) 1415 ±81 1386+^14 

0-074!;}"|^ 0.042 ±0.006 

(F)(10^ergs-'cm-^) 9.05!ff 8.33![^^ 



Notes. Quantities and definitions are as in Table |9] which gives our adopted values, determined assuming circular orbits. Here we do not list 
parameters that are effectively independent of the eccentricity. 

son et al. (2012), and adopted circular orbits, after a significance 
test on the eccentricity measurement. We note that observing the 
secondary transit detection is a very efficient way to characterize 
the shape of the orbit, that would need several tens of RV values 
to obtain an equivalent precision (e.g. Boisse et al. 2009, Hus- 
noo et al. 2011). Another clue to constrain a small eccentricity 
is the measurement of the RVs of planetary lines during the tran- 
sit, but would currently need brighter hosts to be significantly 
detected (Snellen et al. 2010, Montalto et al. 201 1). 

The Rossiter-McLaughlin effect measured via spectroscopy 
during the transit allows the measurement of the sky-plane pro- 
jected angle between the stellar spin-axis and the planetary or- 
bital momentum vector. With the hypothesis that the projected 
angle is close to zero, the expected semi-amplitude of the effect 
is of 26 and 34 m s ' for HAT-P-42 and HAT-P-43, respectively. 

We note that the transit of HAT-P-42b is far from equatorial 
(b = 0.577^y jg5), which facilitates the measurement of the pro- 
jected angle, since there is a smaller correlation between the stel- 
lar u sin / and the projected spin-orbit angle. According to Winn 
et al. (2010a), planets around stars with Teff* < 6250 K are pref- 
erentially aligned (except the systems with longest timescales for 
obliquity damping, e.g. HD80606). Additionally, the (old) age 
of the host stars as well as the (low) mass of the planets also in- 
crease the chances that both planetary systems are well-aligned 
(Triaud 201 1, Hebrard et al. 201 1). We appHed Eq.2 of Albrech 
tet al. (2012) where they estimated the tidal timescale to align a 
planet with the stellar equator from a calibration of binary stud- 
ies. These results also favor alignment for both HAT-P-42 and 



HAT-P-43 (respectively tqe — 247 and 558 to be reported in 
their Fig. 24). 

To conclude, hot gas giant planets are surprisingly diverse, 
and the physical reasons behind the diversity and the observed 
correlations are not well understood. It is likely that the physical 
properties of such planets depend on a large number of parame- 
ters. Hence, there is a continued need to broaden the sample of 
well-characterized transiting planets in order to develop a com- 
prehensive view of exoplanet formation, structure and evolution. 
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